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INTRODUCTION 


OVERVIEW 

This grant funded the development of the thermal characterization 
and analysis of two Space Shuttle power subsystem components: (1) a 
proposed APU line heater, and (2) the fuel cell coolant loop. The 
principal investigator was Dr. J.H. Morehouse who conducted the 
majority of the work on-site at Johnson Space Center (JSC) with the 
Power Branch (EP5) during July and August 1988. 

BACKGROUND 

The hiatus in the Space Shuttle (Orbiter) program provided time 
for an in-depth examination of all the subsystems and their past 
performance. Specifically, problems with reliability and/or operating 
limits were and continue to be of major engineering concern. 

The Orbiter Auxiliary Power Unit (APU) currently operates with 
electric resistance line heaters which are controlled with thermostats. A 
design option simplification of this heater subsystem is being considered 
which would use self-regulating heaters. A determination of the 
properties and thermal operating characteristics of these self-regulating 
heaters was needed. 

The Orbiter fuel cells are cooled with a Freon loop. During a 
loss of external heat exchanger coolant flow, the single pump circulating 
the Freon is to be left running. It was unknown what temperature and 
flow rate transient conditions of the Freon would provide the required 
fuel cell cooling and for how long. 


OBJECTIVE 

The overall objective of the proposed work was the development 
of the thermal characterization and subsequent analysis of both the 
proposed self-regulating APU heater and the fuel cell coolant loop 
subsystem. The specific objective of the APU subsystem effort was to 
determine the feasibility of replacing the current heater and thermostat 
arrangement with a self-regulating heater. The specific objective of the 
fuel cell coolant subsystem work was to determine the transient coolant 
temperature and associated flow rates during a loss-of-extemal heat 
exchanger flow. 
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PART I: SELF-REGULATING HEATER 

The orbiter Auxiliary Power Unit (APU) currently uses electric 
resistance heaters controlled by thermostats to heat the fuel and exhaust 
lines. Keeping these lines at an elevated temperature is necessary to 
prevent freezing and subsequent blockage of the lines, especially when 
the APU is on stand-by during a space operation and the system is 
exposed to space vacuum and low temperature conditions. 

For reliability, the line heaters and thermostats are "doubied-up" 
for each heated line, thus increasing the system complexity. It is also 
necessary to have a large sensing subsystem operating to ensure the 
heaters and thermostats are all operating. In order to reduce the 
complexity of both the heater-thermostat subsystem and the sensing 
subsystem, it was proposed to investigate resistance heaters which are 
self-regulating (self-limiting) to maintain a given temperature. 

BACKGROUND THEORY 

Many materials exhibit the feature of increasing electrical 
resistance with increasing temperature; however, the change in resistance 
is usually fairly small for small (tens of degrees) temperature changes. 
RAYCHEM Corporation has developed an irradiated polymer material 
which exhibits order of magnitude electrical resistivity increases with a 
10° C increase (in the range of room temperature). 

This RAYCHEM material has been proposed as a self-regulating 
or self-limiting electrical resistance heater. The heater is composed of 
two electrical ly-conductive "bus bar" wires on either side of the self- 
limiting material (irradiated polymer), with a voltage difference applied 
across the two wires. As shown in Figure 1, the wires and self-limiting 
material are encased in an insulating material. 

The self-limiting resistance heater works on the same principle as 
other resistance heaters. That is, the heating power is given by; 

Heating Power = I 2 R, 

where I is the current and R the electrical resistance. As shown in 
Figure 2, a voltage (V) is applied across the "bus bar" wires, and the 
current that flows through the self-limiting material is determined by the 
material’s resistance (which is a function of the material’s temperature) 
as determined by the relation: 

V = IR. (2) 

By substituting, the heating power can be expressed as a function 
of the applied voltage and the material’s resistance: 

Heating Power = V*/R ( 3 ) 

Thus, with the resistance directly related to temperature, it is seen that 
the heating power will vary inversely with the temperature of the heater: 
high power at low temperatures, low power at high temperatures. 
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Figure 1 . Schematic of RAYCHEM Corporation self-limiting electrical 
resistance heater. 
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Figure 2. The equations involved in the heating power illustrating the 
inverse rationship between resistance (temperature) and 
heating power. 



The material appears "programmed" with a built-in rheostat 
which causes the material to behave as if it were thermostatically 
controlled. If an environment existed which caused the heater’s 
temperature to decrease, the heater power would increase, which would 
cause the heater temperature to increase. The heater temperature would 
finally stabilize at the point where the heater power equalled the heat 
loss from the heater. The key to this self-regulating heater is that the 
resistance, and thus heater operating range, varies from very low to very 
high values over a very small temperature range. Thus, the heater 
maintains its own temperature within a small range regardless of 
environmental temperature; this is the "thermostat-less" control desired. 

TEST RESULTS 

As discussed, the key to this self-limiting heater concept is a large 
resistance change over a small temperature range. A sample of the selfr 
regulating heater was obtained and a series of tests were run with the 
following objectives: 

• Demonstrate the self-limiting feature of the heater to reach 
equilibrium (room temperature range); 

• Determine the heater power generated per unit length as a 
function of heater temperature; and 

Determine the material’s electrical resistivity as a function of 
the material’s temperature. 

The experimental set-up is shown in Figure 3 where the heater 
was placed inside the environmental chamber taped to an aluminum 
block which acted as a thermal source/sink for the heater. 

Thermocouples were attached external to the encasement material and 
thus did not actually measure the temperature of the temperature- 
resistance material itself. A DC power supply applied a constant 27 
volts to the heater (RAYCHEM instructions) and the amps were recorded 
when the heater was energized. 

Self-limiting Equilibrium Test 

The first test performed was to see if the heater wire would 
actually "self-limit". The test consisted of placing the heater wire 
(attached to the aluminum block) into the environmental chamber which 
was set at 32° F (0° C). After the heater wire and block were chilled to 
32° F and in equilibrium with the chamber, the heater was energized. 

As is shown in Figure 4, the temperature of the heater rose with 
time and reached a "self-limit" of 37.5° F after 5 minutes. The power of 
the heater at this condition was 10W per foot of length. Thus, the 
heater demonstrated the ability for "self-limiting heating". 

Heater Power versus Temperature Test 

The above test at one specific temperature provided the power 
output of the heater for a given heater temperature (10W at 37.5° F). 
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A schematic of the experimental test set-up. The BEMCO 
environmental chamber is located in BLDG. 352. 
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For engineering design purposes with this particular heater, it is 
necessary to know the heater power over a range of heater temperatures. 

By performing the above described equilibrium test with a range 
of environmental chamber temperatures, the heater power performance 
as a function of heater temperature was determined by measuring the 
equilibrium current for each temperature. Figure 5 presents this heater 
performance information. It should be noted that the power 
performance curve shown is specific to the particular heater wire 
(geometry, dimensions, encapsulation material, etc.) tested. 

Electrical Resistivity Determination 

In order to use the power versus temperature information derived 
from the above test, it is necessary to derive a more general, non-design 
specific material property. Electrical resistivity is the material property 
which will allow the design of other heater geometries with other power 
levels for a given temperature. 

The resistivity of the self-limiting material was calculated from 
the voltage and current measurements made when the heater reached 
equilibrium at a given temperature. The resistance was calculated using 
Equation 2, and then the resistivity was determined using the equations 
in Figure 6 and the heater dimensions. The result of these calculations 
at various temperatures is presented in Figure 7. The resistivity variation 
over a 50° F (28° C) range is seen to be two orders of magnitude. 

CONCLUSIONS 

The tests and calculations demonstrate that the RAYCHEM 
material has the electrical resistivity temperature dependence necessary to 
act as a self-limiting (self-thermostated) resistance heater. Use of this 
heater could lead to significant reductions in heater and associated 
sensing subsystems complexity. 

The determination of the resistivity curve allows the physical 
design of heaters with power characteristics determined from thermal loss 
calculations for specific applications. It is recommended that the 
operating characteristics of this self-limiting heater be examined when 
applied to a transient environmental and internal load, such as an 
insulated APU fuel line during start-stop operations. 
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Figure 6. Electrical resistivity relationship to material electrical 
resistance and heater geometry. 
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PART 2: FUEL CELL COOLANT SYSTEM 

The Orbiter fuel cell powerplants generate heat and operate at the 
proper temperature via cooling with a circulated coolant (FC40). This 
circulating coolant rejects the fuel cell heat to the Orbiter cooling system 
through an external heat exchanger. The question has been posed as to 
how long the fuel cell power plant could operate if the Orbiter cooling 
system failed and could not be used as the heat sink. With the loss of 
the heat sink, the fuel cell would ultimately reach a high enough 
temperature that it would be unsafe to continue operating it. However, 
it was unknown how long it would be possible to operate the fuel cell 
before reaching an unsafe operating temperature. 


FUEL CELL TESTS 

Testing for the loss-of-coolant accident profile with a fuel cell 
powerplant was performed at NAS A- Johnson in 1987-88. The tests were 
run at TTA using a prototype 2-stack powerplant, not the Orbiter 3-stack 
system. Additionally, the TTA tests did not try to exactly match the 
external heat exchanger/Orbiter cooling system layout as is found on the 
Orbiter. 


A series of tests were run where the loss-of-cooling temperature 
transients were examined for various operating power levels, both before 
and after the loss-of-coolant. Several of the fuel cell coolant system 
temperatures were monitored during the accident transients, with the fuel 
cell being operated until the temperature of the coolant leaving the stack 
(STKOUT) reached 250° F. Appendix A contains tabular and graphical 
data from some of these 2-stack fuel cell tests. 


PROBLEM AND APPROACH 

Even after the testing of the 2-stack fuel cell powerplants, it was 
still unknown how long an Orbiter 3-stack fuel cell powerplant could be 
operated safely without the heat sink. It also was not possible to test a 
3 -stack system as a "spare" powerplant did not exist. Thus, some 
method had to be developed to relate the 2-stack test data to 3 -stack 
system response to loss-of-cooling. 

The approach used was to develop a computer simulation and 
simulate the accident transient performance characteristics of the Orbiter 
fuel cell powerplant coolant system. The computer model of the coolant 
system was to be verified/validated by the following: 

1) Comparing steady-state simulation results for a 2-stack fuel 
cell with the manufacturer’s (IFC) 2-stack computer model 
results; 

2) Comparing to the TTA transient test results on the 2-stack 
system (Appendix A); and 

3) Comparing to the Orbiter 3-stack system steady-state values. 
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Only after the computer model successfully simulated the above three 
situations would be simulation of the 3-stack Orbiter system transient 
performance to performed the reported. 

THE MODEL 

The computer model of the coolant system involved almost all 
major subsystems of the fuel cell power plant. The schematic of the 
Orbiter fuel cell power plant is shown as Figure 8, with the coolant 
system on the right hand side of the figure. It can be seen that the 
coolant system is directly linked thermally to five major components: (1) 
the fuel cell, (2) the H 2 and 0 2 preheater heat exchangers, (3) the 
external heat exchanger (the heat sink), (4) the condenser for the H 2 and 
H 2 0 mix leaving the fuel cell, and (5) the startup/sustaining heaters (not 
used during the loss-of-coolant accident). 

The flow rates of H 2 , 0 2 , and coolant, the masses and heat 
capacities of all components, and the thermostated and/or pressure 
controlled operating modes, all had to be found, both for the Orbiter 
3 -stack and 2-stack test systems. Much of the data can be found in the 
manual United Technologies FCR-0216, dated 9 April 1976. 

The computer model is named FCPSIM (fuel cell power plant 
simulation) and a complete listing is given in Appendix B. The code is 
fairly well documented using internal comments, and the various 
temperatures, rates of heat transfer, flow rates, and other operating 
conditions are listed and defined. The computer code has been fitted 
with an interactive input format with nine parameters to be set (or 
default value used): 

1) Fuel cell power level (in watts) 

2) Fuel cell type (2-stack or 3-stack) 

3) 0 2 supply flow rate (pounds per hour) 

4) 0 2 supply temperature (in 0 F) 

5) H 2 supply flow rate (pounds per hour) 

6) H 2 supply temperature (in 0 F) 

7) Environmental temperature (in 0 F) 

8) Environmental condition (vacuum or air) 

9) External heat exchanger temperature (°F) 

Figure 9 is a listing of this interactive input format with default values 
shown in parentheses, the choice of fuel cell power of 4100 watts, the 
simulation system time step chosen as 10 seconds, and the printout time 
step chosen as 300 seconds. The output data file (FCPSIM. OUT) 
format is presented in Figure 10 and shows the various system 
parameters which are calculated by the program. The output data 
acronyms are listed and defined within the code, as previously 
mentioned. 

SIMULATION RESULTS 

As described earlier, the verification of FCPSIM involved 
comparisons with manufacturer’s code results for both 2- and 3-stack 
system, and 2-stack transient test results. The 2-stack comparisons, 
steady-state and transient, were able to be done, but the 3 -stack 
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Orbiter Fuel Cell Powerplant 



Figure 8. Schematic of the Orbiter fuel cell power plant with the 
coolant loop on the right-hand side. 
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Figure 10. Example listing of FCPSIM.OUT data listing for two times 
(5 minutes and 10 minutes). 



comparison was not done since manufacturer’s data could not be 
obtained. The results of simulations with the 2-stack and 3-stack system 
parameters is presented below. 

Steady-state 2-stack Comparison 

The initial comparison was between the FCPSIM code and the 
manufacturer’s code (IFC) for steady-state operation of the 2-stack fuel 
cell power plant. Figure II presents the side-by-side comparisons for 
various system parameters, both for low power (2 kw) and high power 
(12 kw) operation. Very good agreement is seen in almost all areas 
except for the coolant flow through the condenser. This coolant flow 
"problem" in the condenser was never satisfactorily understood or 
reconciled. 

Transient 2-stack Comparison 

Figures 12, 13, and 14 present the comparisons between the 
transient temperatures predicted by FCPSIM code and those measured in 
the TTA tests and a 2-stack fuel cell power plant operating at 4. 1 kw, 3 
kw, and 2.5 kw, respectively. Good agreement was produced except at 
the lowest power levels (2.5 kw and below). The major problem appears 
to he in the extremely small amount of excess heat generated by the fuel 
cell, which causes large variations in the temperature changes over time. 
Alternately, the thermal capacitance of the external heat exchanger may 
be chosen as being larger than it "actually" appears. 

The initial (5 to 10 minutes) variations in FCPSIM temperatures 
STKIN and TCE arise from the previous problem of 
controlling/providing the proper flow through the condenser. 

Purge Flow 


It is possible to send a large H 2 flow at low temperature through 
the fuel cell. This low temperature flow is exhausted (purged) overboard 
and acts as a supplementary fuel cell coolant system. Figure 15 shows 
the significant effect this purge flow can have on extending the operating 
time for a given fuel cell power level, from less than 30 minutes to over 
45 minutes with purge. 

CONCLUSIONS 

The FCPSIM code simulating the fuel cell power plant coolant 
system appears to work well at fuel cell power levels above 3 kw, but 
below this level the code does not agree with transient test data well. 
Further investigation is warranted since the lower power levels would be 
expected to be used if an accident involving the coolant sink did occur. 
Also, the 3-stack steady-state FCPSIM code to IFC code still remains to 
be done. 

Of particular note is the fact that using "purge" flow rates with 
the fuel flows would significantly extend the operational time of the fuel 
cell during a coolant loss transient. 
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Figure 12. Comparison of transient temperatures from FCPSIM code 
and TTA test for 2-stack system at 4. 1 kw with loss of 
external cooling sink. 



TEMPERATURE CF) 


HEAT EXCHANGER FAILURE 
RUN*5. 3KW. 62 MIN 


STKOUT 



♦o® TTA TCVT" 


COOLING 
RE^ 


5 2 


TIMf (min) 


Figure 13. Comparison of transient temperatures from FCPSIM code 

and TTA test for 2-stack system at 3 kw with loss of external 
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Figure 15. Comparison of transient stack-out temperatures for 3-stack 
system with and without purge during loss of external 
cooling sink at 4.1 kw (using FCPS1M code). 
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APPENDIX A: 

NASA FUEL CELL TEST DATA 



X-708 DATA 



TIME 

(M«N.) 

0 

1 

2 

3 

4 


S TK IN 
(DEG. F) 

178.70 

178.20 

177.90 

17 8.20 

178.20 


177.90 


177.90 


178.70 

178.40 

178.70 


178.90 


179.10 


C 


2.5 KW 


STK OUT 
(DEG. F ) 

195.00 

194.0 0 
193.30 
192.80 
192.60 


192.60 


192.30 


192. 30 

192.60 

192.80 


193.10 

193.10 


TCE 

( DEG. F) 

146.20 

147.20 

148.40 

148.40 

148.40 


148.60 


149.90 


1 49.90 

150.60 

151.10 


151.80 




162.50 


164.30 

t — — 


165.70 


180.10 

199.40 

167.40 

180.90 

199.20 

168.70 



182.60 


185.30 


187.00 


188.70 


190.10 


191.40 


193.30 


194.50 

196.00 

197.00 


198.40 
199.20 

200.40 


201.60 


202.60 


203.60 


204.50 


2 05.80 


207.00 


198.90 


199.20 


199.20 


200.10 


200.90 


201.40 


202.10 


202.80 

2 03.60 

204.50 


205.50 

206.50 

207.50 


208.70 


209.70 


210.90 


211.60 


212.80 


213.60 



174.80 


179.10 


182.30 


184.30 


186.20 


187.70 


189.90 


191.4 0 

1 93.30 

194.30 


195.80 

197.00 

197.90 


199.20 


200.40 


201.60 


202.60 


203.80 


204.80 
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c 


2.5 KW 


STK OUT 


229.90 

230.70 

231.40 

2 32.40 

233.10 
233.80 


234.60 

235.30 

236.00 


236.80 

237.70 


238.20 

239.00 

239.70 

2 40.70 
241-20. 
242.60 

243.10 

183.50 

177.90 


223.60 

224.30 


225.50 

226.50 


227.00 


228.20 


229.00 

220.40 

229.90 

221.60 

230.70 

222.40 

231.60 

223.30 


232.40 



233.30 

225.30 



234.30 

226.30 



235.10 

226.80 



235.80 

227.70 




237.70 

238.40 
239.0 0 
239.90 

240.40 


241.40 


242.10 


242.60 


243.60 


244.30 


245.10 

245.6 0 

246.50 

247.50 

248 .00 

249.50 

250.00 

216.70 

207.00 


229.20 

230.20 
230.70 
231.90 
232.40 


233.30 


234.10 


234.80 


235.50 


236.30 


237.00 

237.70 

238.50 

239.20 

239.9 0 

2 41.40 

241.90 

165.20 

155.70 




































































X-708 DATA 


STKOUT 





179.60 


177.90 

1 79.10 

179.40 

179.10 

179.40 


179.60 
1 79.90 
180.10 
180.40 


180.60 

180.90 


180.90 


180.90 


1 89.20 

191.10 

193.10 
195.00 
196.70 


198.40 

199.90 


201.40 


202.80 


204.30 


205.50 


207.00 


208.20 


209.90 

210.90 

212.60 

21 3.6 0 
215.00 

216.30 


DEG. F 

) i 

( 

DEG. F 


197.20 


196.50 

196.00 

196.00 

196.20 

196.50 


196.70 

197.00 



207.00 

208.20 


209.70 


210.90 


212.80 


214.10 


215.50 


216.70 


218.20 


220.70 

221.90 

223.30 

224.30 


146.20 


148.10 

1 50.10 

149.10 
150.30 


150.80 


151.60 

152.00 



197.90 

154.70 

197.90 

155.70 

198.70 

157.90 

199.40 

160.10 

200.40 

162.10 

200.10 

164.00 

199.70 

166.90 

199.40 

168.40 

200.40 

170.10 


195.00 


197.00 


198.70 


200.10 


201.60 


202.80 


204.50 


205.80 


207.50 


210.40 

21 1 .60 

213.30 

214.30 




175.20 


180.60 


185.00 


187.20 



220.20 

228.50 

221.10 

229.40 


218.50 


219.40 






41 

218.70 

227.00 216.70 
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X-708 DATA 



8 

0_ 

8 

1 

8 

CM 

8 

3 

8 

4 

8 

5 

8 

6 

8 

7 

8 

8 

8 

9 
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X-708 DATA 



TIME 

(MIN.) 

0 

1 

2 

3 

4 


STK IN 
JDEG.F1 

17 8.70 

177.70 
_ 179.40 

178.90 
180 JO 

180.40 

180,60 

179.90 

180.60 

180.6 0 

18 1.1 0 

18 3.10 

189.60 

193.10 

196.00 

198.40 

201.10 

203.30 

2 06.00 

208.20 

210.60 

212.60 

215.00 

217.00 

219.40 

2 21.40 

22 3.30 

225.80 

227.50 

229.70 

231 .60 

233.80 


4.1 KW 
STK OUT 
(DEG. F) 

1 99.40 

190.20 

197.00 

197.20 

198.40 

199.9 0 

201.60 

201.90 

201.60 

202.10 

203.1 0 

203.80 

204.80 

206.20 

208.40 

210.40 

212.80 

214.80 

217.50 

219.40 

222.10 

223.80 

226.50 

228.20 

230.40 

232.40 

234.30 

236.80 

238.20 

240.70 
242J0 

244.60 


TCE 

(DEG. F) 

148.90 

148.10 

150.30 

152.30 

1 56.20 

156.70 

1 59.60 

161.80 

167.20 

170.80 

174.30 

175.70 

183.10 

188.70 

192.10 

194.50 

197.70 

200.10 

203.10 

205.30 

208.00 

210.20 

212.60 

214.50 

217.20 

218.9 0 

220.90 

223.30 

225.30 

227.50 

229.2 0 

231.60 


237.50 

239.20 

218.90 

199.90 




248.50 

250.00 

2 45. 60 

234.30 




235.30 

237.0 0 

205.00 
183.10 
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Dr. Jeffery Morehouse 
Department of Mechanical Engineering 
University of South Carolina 
Columbia, S.C. 29208 

Dear Jeff, 

I could not find any of the X-708 test data or records of the 
previous work on the fuel cell simulation program. All that 
information was stored in my old office over the past year. Nobody 
appears to know what was done to the records when the new people 
in EP5 moved into that office. I reconstructed the fuel cell stack in, 
stack out, and condenser outlet temperatures from some plots stored 
in the computer. I also included the default input parameters for the 
simulation program. I know that some of the flow rates and heat 
exchanger temperature values will be off. During the X-708 test 
program the fuel cell was purged during the tests for failure of the 
fuel cell heat exchanger. You will just have to guess what those flow 
rates were. If you have any questions call me at (713) 483-9048. 

Sincerely, 

- tr Jr 

Howard A. Wagner 



APPENDIX B: 

FCPSIM COMPUTER CODE LISTING 



& kkk PROGRAM FCPSIM (FUEL CELL 
C A A A I N p ij x S 


S I M U L A T 1 0 N ) 




10 

15 


20 


30 


: x a: 


: ; ; p 


; R E 


' 5 l * EC POWER •: 1 y 'v '•} , 

: 2 U P P L Y F L 0 w R a t E : 

I nDa-uRs '?5or«> vsx^ss ° enmirom2rntai.' • 

* ; Cy A Y; x lY / 5 : Y' 9 C EXTEftNA; - heat exchange® 

( :-'0 . 0 > ;/",/ , HIT RETURN WHEN FIN]'-‘HPn rn 

■ EXE C U X E P R 0 G R A M ' / ~ 

shad <5,25) NF '... AG 
I F ( N E L A G . E S , 1 ) T H E N 
W R I X E ( 6 , 1 0 ) 

FORMAT-; IX, 'ENTER FUEL CELL POWER IN WATTS"- 
READ (5, 155 PF C 
FORMAT (F8. 3 > 

GOTO L 
END IF 

IF < NFL AG . BQ.2 ) THEN 
WR I T E < 6 , 2 0 ) 

FORMAK IX, 'ENTER FUEL CELL TYPE'/5X ' 

& '2: 3 SUBSTACK' ) 

R E A D ( 5 , 2 5 ) N E W 
FORMAT (ID 
G 0 TO 1 

e n d f 


2 S j B S T A C K ■' , 10 X 


IF ( NE W , SQ . 1 ) OLD*. TRUE. 

Dr (NEW. EQ.2)0LD=. FALSE. 

IF (NFLAG.EG.3)THEN 
WRITE <6, 30) 

FORMAT (IX, 'ENTER 02 SUPPLY (REACTANT P 
i 'LB/HR') 

READ<5, 15) W02 
GOTO 1 
END IF 

I F ( N F L A G . E Q . 4 ) T H E N 
W R I T E ( G , 2 5 ) 

FORMAT (IX RENTER 02 SUPPLY TEMPERATURE 
READ (5,15) T02IN 
GOTO 1 
END IF 

IF (NFLAG. ECUS) THEN 


U S 


IN 


PURGE) FLOW RATE 


) 


OF POOR 


O A («*e“ 

2 if •/. Cc i_ 





WRITE; G , 40) 

4.0 FORMAT* I X, 'ENTER H2 SUPPLY (REACTANT PLUS PURGE) RATE I 

R E A D ( 5 , 1 5 ) W H 2 

E N D I E 




b N .U i. 

IF<NFLAG.Ea.9>THEN 
w R I T £ ( 6 , 6 0 ) 

F 0 R M A T < 1 X » ' E N TER E X T E P N A L H E A T E X C H A N G E R TEMPE R A T U R E 
R E A D < 5 , 1 5 ) T I" A X 
G U T 0 1 

W K 1 T E ( 6 , 6 5 ) 

F 0 M A I (. I X , 'ENTER C A L C U L A T 1 0 N TIM E S T E P 3 I Z E I N 3 E C 0 N D S 
R H A D ■ 5 , 15 ) D TIME 
B 1 1 M E D T I M E , ’ 3 6 0 0 . 0 
WR ITE < 6,70) 

FORMAT ( IX, 'ENTER PRINTOUT TIME STEP SIZE IN SECONDS' > 
READ <5, 15) PI IME 
?TIME=PTI ri E / 360 0.0 


■■ E C 


C kkk OUTPUT HEADINGS 
C 

OPEN ( UN IT = 3 , F ILE- ' FCPSIM. OUT ' , STATUS- ' UNKNOWN ' ) 

I F ( 0 L D . E Q . .TRUE. )THEN 
STACK- '2 SUBSTACK' 

ELSE 

STACK* '3 SUBSTACK' 

E N D IF 

L F ( Z . L E . 0 . 5 ) T H E N 
VAC- ' VACUUM ' 

ELSE 

VAC- 'AIR 
END IF 

WR ITE ( 3 , 90 ) RFC , STACK , VAC , W02 , T02 IN , WH2 , TH2 IN , TENV , TEHX 
90 FORMAT ( IX , ' FUEL CELL WATTS: ' , F 7 . 1 , 3 X , A 1 0 , 3 X , A 6 / I X , 0 2 SUPPLY 
i 'FLOW RATE (LB/HR): ' ? F5.2,3X, 02 SUPPLY TEMPERATURE (DEG F) 

X F 6 . 1 / 1 X , ' H 2 SUPPLY FLOW RATE (LB/HR): ' , E5 . 3 , 3X , ' H2 SUPPLY 

X 'TEMPERATURE (DEG F :• : ' , F 6 . 1 / IX, 'ENVIRONMENTAL TEMPERATURE 

■i '(DEG F): EG. 1/ IX, 'EXTERNAL HEAT EXCHANGER TEMPER ' y 

4 'ATURE (DEG F): ' ,F3. l//> 

WRITERS, 100) 



1 00 F G R M A T ( 2 X , 'TIME T S K 0 

X '■ QFCRH CURR UFC 

t ' T H 2 0 U T T C P H H Q p H H 2 

& '' WCB1 WCB2 WCEHX 

& ■' TCE QC VEN 2CC' 

% ' T H I hi !••: 2 I T H " •*; T 


irC _ QQSN QCFC QFCHN ’ 

T C P H 0 ' / 2 X , ' T 0 2 0 U T QPH02 n P M fj f 
q p -i r f-J T C P Q c P E N w c ; / o X 

* •- 2 T 2 H > Q H X 2 N Q £ h ;< 

2 X ' T C C 0 Q C :j H D 1 c P £ N 


P H 2 H 


A A A Ui A X E .E 


H E A X A N D L A X E N 


H E A T L 1 0 0 6 „ 4 

*AA thermal properties 

S I G ~ 0 , 1714 E -• 8 
EM I S •“ 0 . S':.: 

H c u n y i . a 


"! E H i. 


AAA THERMAL MASS OF C 0 M P 0 N E N T S 
AAA FUEL CELL 


IF (CL D, 

EQ. 

.TRUE. ) C M E L “ 2 6 . 5 A 0 . 7 

IF COLD. 

EQ . 

. FALSE . ) CMEC~39 . 75 AO 

OXYGEN 

PREHEATER 

CMPHO-O 

. 16 


A P H 0 = 0 . 

23 


E P H 0 2 = 0 

. 9 



Aa.A HYDROGEN PREHEATER 

C M P H H = 0 . 10 
A P H H - 0 .146 
HPHH2-0 . 9 

AAA COOLANT PUMP 

WKCP-41Q . 0 
0 M C P ~ 0 . 42 
A C P = 0 . 795 

A*A EXTERNAL HEAT EXCHANGER 

ifcqld.eq. .true. )then: 
CHEHX=7.09A0.5 
A E H X = 6 . 03 AO 3 


Op FQOR. Q*u 



ELSE 

CMEHX-2 . 19 
A E H X = 0.904 


!'. A- 4- tV 


NT POL 


, : / :• -j s x a i n : g heater s 



C A k k I N 1 1" I A L. V A L U E S A A A A A A A A A A A A A A A 

C A a A F U E 1. C E L L, 

C 

TSTQP= 11 0.0/60,0 

I 3 T A R T = 2 0 . 0 / 6 0 . 0 

T 3 K I N = 1 8 0 . 7 

T F C = 0 , 0 0 4 7 A P F C + 1 7 8 . 6 

TSK0-IFC-2.Q 

UC = 4 25 . 0 

V F C = 3 0 . 0 

W 0 2 F = 0 , 6 7 A P F C / 1 0 0 0 . 0 
W H 2 F = 0 . 0 8 5 A P F C / 1000.0 
WH2I-WH2 
WQ2I=W02 
C 

C AAA 02 PREHEATER 
C 

T 0 P H 0 = T S K 0 - 2 . 0 
Q P H 0 2 -• 1 0 0 . 0 


C AAA H2 PREHEATER 
C 

T C F H H = T C P H 0 - 2 . 0 

C AAA COOLANT PUMP 

X C P = T C ? H 0 + 1 . 0 
C 

C AAA CONDENSER CONTROL 


C 


T C E = 1 5 1 . 7 
T C E P P = I C E 
TCEP-TCE 
WCB2=0 . 02APFC 

CONDENSER 

ID P=0. 00 3 APEC +154 
TCCQ~ .1 77 . 7 


VALUE 


0 


AAA 


TCCOP-T CCG 
TCC 1 = 1 13.0 

iccip=rcc i 

'4 H 2 0 R ~ 6 .. 15 
N H 2 S E C 1 0 „ 2 



CURR-PFC/VFC 
IF ( OLD . EQ . . TRUE . ) THEN 
I FIGURE. LE. 131 ) THEN 
S L 0 P E = - 0 „ 0 1 7 1 
Y I N T = 3 2 .67 
ELSE 

S L 0 ? E ~ - 0 . 01093 
Y INT=31 . 87 
END IF 


- — 




\ 



o o 


ELSE 

IF(C‘JSR.LE. ISO 5 THE : 
S L OPE = - 0 . 0 1 0 4 


E N ft T F 



C 

C 

c 



3 E 


: h L- 


W02ACP02H ) 


Q G 2 N “ E C K W A < 1 . 0 

WH2F-CURRA2 . 35 
WH2QF~CURRA2 . 3 
'I* 0 2 0 U I - T 0 2 I N + Q 
IJH23 = WH2QR + WH2GF 
W H 2 T -- W H 2 F + W H 2 F E C 
T H H P = I C E + 2 „ 0 
! F C P -- T F C 

U F C K E — ~ < WH2QRA L P H 2 0 + W H 2 A C P H 2 ) A T H HP~W02 A C P 0 2 A X 0 2 0 IJ I + ( W 
& + W H 2 R E C A C P H 2 ) A T E C P + ( W 0 2 2 02 F ) A C P 0 2 A T F C P 

Q F C E N = 1 . 163E-10A < ( TSK0+460 . 0 ) AA4- ( TENV+460 
% ( 0 . 2 4 9 A A L 0 G .1 0(2,1 b E 6 A C TSKQ-TENV ) ) - 0 . 27 ) A 

TFCP- ( S3 GEN- Q F C R E - Q E C E N + W C A C P C A < 2 + TSK IN) / < W 
T F C P - ( T E C P + T F C ) / 2 . 0 

T F C ~ ( G G E N - Q F C R E - Q F C E N + W C A C P C A < 2 . 0 + T S K I N ) + C 
% / ( U C A C ? C + C M F C / D T I M F. ) 

DELTFC - ( QGEN-QFCRE-QFCEN + ECACPCA < TSK IN-TSK 
T F 0 = T F C P + D E L T F C 


, 0 ) A A 4 
(TSKO- 
ICACPC) 


) + Z A 0 . 

TENV) 


1 M F C A T F C P / D T I 
;0) ) A (DTI ME/ C 


H 2 0 A C P H 2 ! 
0 0 3. 3 6 A I 0 ; 

H E ) 

M F C ) 


T D R 0 P = ( 4 . OAWC/1 100 . 0 ) + ( 4 . 0APFC/120C0 , 0 ) 
TS K 0 = T F C - T D R 0 P 
QCFC=WCACPCA( TSKO-TSKIN > 


C AAA 02 PREHEATER 

C 

AAA PARAMETERS 


C 


0 

C 

r 

c 

c 


CMPHO (MC) OF 02 PREHEATER 

APHQ SURFACE AREA (0.28 SO FI) 

WD2 FLOW RATE OF 02 (LB/HR) 

TENV ENVIRONMENTAL TEMPERATURE (DEG R) 

EM IS EMMISIVITY OF PREHEATER SURFACE (0.35) 

HCONV CONVECTION COEFFICIENT IN AIR (BTu/HR SQ FT DEG R ■ 

Z FLAG (1.0 FOR AIR ENVIRONMENT, 0 FOR VACUUM ENVIRONMENT) 

TQ2IN 02 SUPPLY TEMPERATURE ; DEG R) 

EPH02 EFFECTIVENESS OF 02 PREHEATER (0.90) 

SIG STEPHAN-BOLTZMAN CONSTANT (0. 1 7 1 4 E - 8 BTU/HR SQ FT DEG R A A 4 
CP02H SPECIFIC HEAT OF 02 SUPPLY 


OF PO 



CK 


;/ 





IC?HOR=TC?HQ 
TSKOP'-TFCP” TDROr 

H y 2 — t P h‘ U 2 A W !J 2 * •_ y *J 2 H A . '■> 0 ~ 1 0 ... 

J p :J ~ N - S 1 . G A r, i Vl I A A p !-i i j * ! . o r- ■" *. ) ) ) 1 •- . H N L f 4 b Q . .) ) A k 4 ■ * t • 

p Li n + p . p - j £ v y 


c 




U P H H E N S I G A E M I A A p Pi h T C P h 0 + 4 b 0 . 0 ) A A 4 ~ ( TEN V + 4 6 0 . 0 ■ A A 4 ) + Z A H 

1 A F* H H A ( i C P ' 3 - T E N V 

I c ? H u P ~ ( - r; ? H H E N. » G p hi ri 2 + Ui r; * c P C A X c p hi 0 ) / ( w C A C ? c ) 

T 3 P H H - ( - Q P i-i HEN - : ? S H 2 L U C A C P C A X C p H 0 ~ 3 M P H H A (TCP H 0 / 2 . 0 - ( T C P A H p + 
% T C: p H 0 R ) / 2.0) / D T I M E ) / ( W C A C P C + C H P hi H / < 2 . 0 A D T I M E ) ) 

T ~ 2 3 U T :: = T H 2 IN + G P A i-i 2 / ( u) i-i 2 * C P H 2 H 


n >v 


C A A. A C l CL A NT PUMP 

u A A A p A R A METER S 
C 

C WKCP POWER INPUT TO PUMP MOTOR (120 WATTS) 

C ACP COOLANT PUMP SURFACE AREA 

C CMCP (MO OF COOLANT PUMP 

(.'! U P E N ~ S I G A E M I + A tj p A >; i T C P H H + 4 6 0 . 0 ) A A 4 — ( TENV + 4 G 0 . 0 ) A A 4 ) + ZAH C 0 N V A 
i ACP A ( TCPHH-TENV > 

T C p p - < ~ Q C P E N + W K C ? + W C A C P C A T C P H H ) / ( WCACPC ) 

TCP~ -QCPEN+WKCP+WCACPCATCPHH-CMCPA ( TCPHH/2 . 0- ( TCPP + 

X TCPHHR ) /2 . 0 ) ./HUME ) / ( WCACPC +CMCP/ ( 2 . OADT IME ) ) 

C AAA COOLANT CONTROL VALUE 


c 

AAA PAR 

AMETERS 



! t 

A 1 

LOWER TEMPERATURE 

SETT I NO 

(DEG F 

C 

A2 

UPPER TEMPERATURE 

SETT I NO 

( D E f 3 E 

c 

B 1 

MINIMUM FLOW RATE 

( L B / H R ) 



B 2 

MAXIMUM FLOW RATE 

( L B / H R > 



A 1 ~ 

1 3 6 . 0 




A2 - 

203.0 




IF (OLD. EG. . TRUE. ') TmEN 
61=425.0 

boo : oo.o 

ELSE 


OF 'POOS QUALiTY 





S 1 =8 • j 0 •» 0 
B2= 1 400 . 0 
END IF 







ELS) 


I F ( T M I X P ... L E . 1 8 1 . 0 > U C B 1 •= W C - 1 4 0 0 

I F •' ( 7 M I X P G T . I 8 3. . 0 ) . A N D . < T M I X p . L T „ 185. 0 ) ) W C B 1 ~ ( W C 
X - 1 H I X P ) / 4 „ 0 

•v L C “ W U ~ j» C B 1 

1 1 \ u C 0 N D • G T . QCO N U C ) +> C = ( G C 0 N D A WUC > / Q C 0 N D C 
V F ( W C G . LT . 0 „ 0 ) w C C :::: W C 
i E w U U . i i E . W C .> W C C “ W C 
w c b i " w c - wcc 

C AAA ENERGY BALANCE ON NIXING 0 A L V E G I V E S X C C 0 
C 

TCCO = ( WCATM IXP-WCB 1 ATCP) /WCC 

C 


4 0.0' 


C 

C- A A A CONDENSER (BACKWARDS THROUGH IT!) 

C 

WH20=WH20F+WH20R 
WH2T=ldH2 + WH2REC 

TH I ~ ( W H 2 A C P H 2 A T H 2 Q U T + W H 2 R E C A C P H 2 A T F C + W H 2 0 A C P H 2 0 A T F C ) (WH2T 
X ACPH2+WH20ACPH20) 

GSENH- ( UH2TACPH2+WH20ACPH20 ) A ■ TH I-TDP •' 

QLAT=HEATLAWH20F 

Q3ENL= ( WH2TACPH2+WH2QRACPH2C) * (TDP-TCE ) 

IF < (I IME.GT.TSTART) .AND. < TOE .. GT . TDP :■ ) THEN 
QSENL = 0 . 0 
Q L A I = 0 . 0 
END IE 

QC0ND=G3ENH+GLAT+GSENL 
ICC I P - T C C I 

G C D E N = 3 I G A E M I S A A C 0 N D A ( ( T i i + 4 G 0 , 0 ) A A. 4 ~ l. T E N 0 + 4 6 0 . 0 ) A A 4 ) + Z -I C 0 N 0 A 
% ACONDA ( TH I-TENV ) 

TCC 1= ( WCCACPCATCCO+QCDEN-QCGNB / / ( WCCACPC ) 

C T CC I" < w : C A C. P C ATCCO+ 3 CDEN - 0. C 0 N. D + CM C C N D A (TCC 0 -TCC 0 F-TCCI ) / 

C % ( 2 . OADTIME ) ) / < WCCACPC -CMCGND/ ( 2 . 0 ADI IMS ) ) 

IF < TCC I . L E . T E H X ) T C C I - T E - X 
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W C B 2 - 0 , 0 
!« C E H X ~ W C C 


C TCCI-TEHX 

C END IF 

■J AAAAAAt^AAAAAAA AN 0 X H HR TRY ! 

W C B 2 ™ f -a C V E N + w C C A C ? C A ( T C C I - 1 E H X ) ) / ( CPCACT C P T E H X ) ) 
T C E - 1 5 3.0-5 . OA W C E 2 / W C C 
I F ( W C B 2 . L E . 0 . 0 ) T H E N 
TCCI=TEHX 
WCB2 = 0 . 0 
E N D I F 

IF ( WCB2 . GE . WCC ) WCB2-WCC 
IF ( TCE .LE.TCC 1 + 2 . 0 )TCE-TCC 1 + 2 . 0 
IF ( TCE .LE. 148.0 ) WCB2=WCC 
IF ( TCE .GT.153. 0 ) W C B 2 0 . 0 
WCEHX=WCC -WCB2 
C 
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C QHXEN 
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C TEHXP 
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EXTERNAL HEAT EXCHANGER TEMPERATURE (DEG S) 

EXTERNAL HEAT EXCHANGER COOLANT FLOW RATE (LB/HR) 

STEPHA N - B 0 LTZMAN C 0 N S T A N T 
EMISSIVITY 

ENVIRONMENTAL TEMPERATURE (DEG R ) 

FLAG FOR E N V I R 0 N M E N T A L lON 0 I T 1 0 N S 
CONVECTION COEFF ICIENT 
3 PECIFIC HEAT 0 E T K > C 0 0 L A N T 
COOLANT PUMP TEMPERATURE (DEG ft) 

(MC) OF T E EXTERNAL HEAT EXCHANGER 
TIME STEP 


R QUALITY 
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[ T I IN .. T H 
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t c r 

TEMPERATURE OF THE 

F 
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(DEG ?> 
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FLOW RATE OF WATER 
SPECIFIC HEAT OF W 

AT 

H R 

HE FUEL CELL <L 

QCO t- j 

HEAT TRANSFER IN T 

H E 

0 0 N D E N C 

E R i E T ! : / H ft > 

T f wT 

HEATL 

TEMPERATURE AT I ON 
LATENT HEAT OF CON 

U E 

NS HR EX I 
N 3 A X 1 0 N 

T (DEG TE) 

G U D E N 

E N V IRONMENTA L H E A T 

L 

0 S S <■ 3 T U 

/HR ) 

T C C 0 

TEMPERATURE OF COOL A 

N T E X I T I 
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w c c 

C i'i C 0 N [i 

COOLANT F L 0 W SATE 
CMC) OF THE CQNPEN 

X H 
S £ 

ROUGH CO 
R 

N l.i E N S E R (LB / H R ) 

TOO OP 

PREVIOUS TEMPERATU 

RE 

0 G 0 L A n T 

EXITING CON DEN 

p p TC 

D T I :1 E 

PREVIOUS CONDENSER 
TIME STEP 

E 

X IT H2 TEMPERATURE 


T C COP- TC 0 0 

QC0ND-WH2TACPH2A ( THI-TCE ) +WH20F AHE A IL+WH2G ACPH20A ( T H I -T 
T C ! :J ~ ( - Q C D E N + Q C C N D + W C C A C P C A ICC I ) / ( W C C A C P C ) 

TCCO- ( -QCPEN+aCOND+wCCACPCATCC I+CMCOND A ( ICCQ+TCC IP 
l ™TCC I > / ( 2 . OADT IME ) ) / ( WCCACPC +CMCQND/ < 2 . 0 A D T IME ) ) 

A A A TO GET TCCO < AND THUS TCCI ALSO) TO ’SEEK* LOW OPERA! [NG 
X C C 0 M = T E H X + Q C 0 N B / ( WCCACPO '• 

IF < ( TCCO . GT . TCCOM ) .AND. ( T S E , G T . 152 .0) ) TCCO= < TCCOM + TCCO ) / 
AAA TO PREVENT TCCO EXITING HOTTER THAN SOURCE 
IF ( TCCO. GE.TH 1-2.0 )TCCQ = TH 1-2.0 
AAA TO BALANCE HEAT FLOW FROM M2 TO COOLANT 'SEE STACK INLET 
QCONDC = UCCACPCA ( TCCO-TCC I ) 
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